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Abstract. We give a necessary and sufficient condition for tlie giobal 
' existence of tlie classical solution to the Cauchy problem of the com- 

pressible Euler-Poisson equations with radial symmetry. We introduce 
a new quantity which describes the balance between the initial velocity 
of the flow and the strength of the force governed by Poisson equation. 

Ph 

[ 1. introduction 

We consider the compressible Euler-Poisson equations: 
(1.1) pt + dw{pv) = 0, 

>; (1.2) vt + v -Vv = -XV^, 

g ; (1.3) A<i> = p, 

(N 
(N 



where {t, x) € K_|- x M^. These are the conservation of mass, Newton's second 
_ law, and the Poisson equation defining, say, the electric field in terms of the 

■ charge, respectively. The unknowns are the mass density p = p{t,x), the 

I velocity field v = v{t, x), and the potential $ = ^{t, x). A is a given physical 

' constant. 

In this paper, we assume that the unknowns have radial symmetry and 
^ \ concentrate on the multi-dimensional isotropic model: 

(1.4) r"-V + 5.(r"-V) = 0, 

(1.5) Vt + vdrv + \dr<l> = 0, 

(1.6) a,(r"-ia,$) = r"-V 

for (t, r) G M+ x M+ with initial data 

(1.7) {p,v){0,r) = {po,vo){r), po ^ 0. 

Here, r ^ denotes the distance from the origin. Now, the unknowns are 
p = p{t, r), V = v(t, r), and $ = <I>(t, r). 

The Euler-Poisson equations arise in many physical problems such as 
fluid mechanics, plasma physics, gaseous stars, quantum gravity and semi- 
conductors, etc. There is a large amount of literature available on the global 
behavior of Euler-Poisson and related problem, from local existence in the 
small ff^'-neighborhood of a steady state [19, 21, 11] to global existence of 
weak solution with geometrical symmetry [8]. For the two-carrier types in 
one dimension, see [25]. The relaxation limit for the weak entropy solution, 
consult [22] for isentropic case, and [15] for isothermal case. The global 



existence for some large class of initial data near a steady state is obtained 
by Guo [14] assuming the flow is irrotational. 

For isotropic model, the finite time blowup for three dimensional case 
with the attractive force, pressure, and compactly supported mass density is 
obtained in [20] , and the blowup for the repulsive case in the similar settings 
is deduced in [23]. In [10], the global existence/finite-timc breakdown of the 
strong solution is studied from the view point of critical threshold. They give 
a complete criterion in one-dimensional case without spatial symmetry and 
with spatial symmetry in one and four dimension. A sufficient condition 
for finite-time breakdown without spatial symmetry is obtained in [7, 6], 
and the complete description of the critical threshold phenomenon for the 
two-dimensional restricted Euler-Poisson equations is given in [18]. 

In this paper, applying the method in [10], we discuss the necessary and 
sufficient conditions for the global existence of the solution to the Euler- 
Poisson equations with spatial symmetry (1.4)-(1.7) in multi-dimensional 
case. 



1.1. The Euler-Poisson equations and the Schrodinger-Poisson sys- 
tem. The Euler-Poisson equations (1.1)-(1.3) are related to the Schrodinger- 
Poisson system via semiclassical limit. The global existence of the solution 
is applicable to the study of the Schrodinger-Poisson system. Consider the 
Cauchy problem 



(SP) 



\u 



h\2 



u"(0,x) = ao{x)e~, 



where h is a positive parameter corresponding to the scaled Planck's con- 
stant. For simplicity, we suppose that Oq and in the initial data are 
independent of the parameter. In the limiting process ^ — ^ 0, following 
WKB type approximation is often considered; 



(1.8) u^{t,x) ~ e*~?^(a(t,x) + hai{t,x) + h^a2{t,x) H ). 

One way to justify this approximation is to employ Madelung's transform 

u^{t,x) = ^Jpf'{t,x)e'~h— 
and consider the quantum Euler-Poisson equations 

dtp'' + div(/V5'') = 0, 

dtVS'' + (VS'^ • V)VS'^ + XVVp = — V 



(QEP) 



avp' 



AK 



(/(0,x),V5'^(0,x)) = (|ao|^V(<^o + /iargao)). 
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0, 



The term (/i^/2)V(Ayp^/yp^) is called quantum pressure. Taking h 
we obtain, at least formally, the Euler-Poisson equations 

dtp + div{pv) = 0, 

dtv + {v ■ V)v + XVVp = 0, 

AVp = p, 

ip{0,x),viO,x)) = i\ao\\V(t>o), 

where p = lim^-^o P^, v = lim/i_»o VS''*. This limit is treated in [26] and the 
convergence of the quadratic quantities 



(EP) 



h\2 



p, hlm{vP-Vu^) ^ pv 



is proved in the sense of Radon measure with some restrictive assumptions. 
Though this convergence suggests that the solution vl^ may have the asymp- 
totics vP' = e^^/^^^/p + o(l)), it is not satisfactory. 

Another way to justify (1.8) is to employ a modified Madelung transform 



(1.9) 

and consider the system 



a e h 



(QEP') 



dtV^^ + {V^^ ■ + XVVp = 0, 



AK 



(a'^(O,x),V*'^(O,x)) = (ao,V0o)- 



It is essential that takes complex value, and so, in general. The 

point is that the system (QEP') can be regarded as a symmetric hyperbolic 
system with semilinear perturbation. It is proven in [2] that this system is 
locally well-posed for ^ /i <C 1 and the solution (a^, "^^) can be expanded 
as 



a + hai + y?a2 + 



Plugging this to (1.9), we obtain WKB type estimate 

(1.10) n'*(i, X) = e'^ (/3o(i, x) + hpi{t, x) + hp2{t, x) + • • • ) 

with f3o = ae**i in a time interval which is small (in general) but independent 

of the parameter. This method is first applied to analytic data ([12]) and to 
Sobolev data ([13]) for certain class of nonlinearities, and it is generalized 
in [1, 9, 16, 4, 17, 2, 5]. We also note that the approximation of the form 
(1.10) leads to some ill-posedness results for the "usual", that is, non-scaled 
nonlinear Schrodinger equations ([3, 24, 5]). One verifies that the principal 
part of the solution (a, ^') solves 



(1.11) 



dta + (V* • V)a + -aA^ = 0, 
dtV^ + (V* • V)V* + XVVp = 0, 

AVp = |a|^ 
(a(0,x),V*(0,x)) = (ao,V.^o)- 
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Hence, we see that p := |ap = |/3oP and v := also solves (EP). 

Either way, the problem of the justification of the global estimate of the 
form (1.10) is closely related to the problem of global existence of the solution 
to (EP). If the solution of (EP) is not global and breaks down in finite time, 
it immediately implies that the WKB type estimate breaks down at the same 
time. The converse is not so clear. In one-dimensional case, this limit and 
the large time WKB type estimate is given in [17] using the result in [10]. 

1.2. Critical thresholds. Before stating our main result, we briefly recall 
the part of the result in [10]. They introduce the notion of critical thresholds 
and give several sufficient conditions for global existence and finite-time 
breakdown in terms of the initial velocity. We restrict our attention to 
the positive time t ^ 0. The necessary and sufficient condition for global 
existence is obtained in the case n = 1 or 4. For a nonnegative integer s, we 
define 

('C([0,oo)) if 5 = 0, 

|(7([0,oo)) n(7^((0,oo)) ifs>0. 

Theorem 1.1 (Critical thresholds in ID case, [10]). Suppose n = 1, A < 0, 
Pq G D*, and vq € D^^^ with vq{0) = for some positive integer s. Then, 
the classical solution to (1.4) -(1.7) is global if and only if 
(1-12) 

vo{R) > -^2\X\rJ^ pois)ds and v'o{R) > -y/2\X\po{R), Vi? > 0, 

where, in both inequalities, we allow the case where the both sides equal zero. 
If po and Vq satisfy (1.12) then the corresponding solution of (1.4)-(1.7) 

satisfies 

p G C2([0, oo), D') n C°^{iO, oo), D'), 

V G C\[0, oo), D'+^) n C°°((0, oo), D'+^), 
$ G C2([0, oo), D'+^) n C°^((0, oo), Z)^+2). 

The solution is unique in C'^{[0,oo), D^) x C^{[0,oo), D^) x C'^{[0,oo), D'^) 
and also solves (1.1) -(1.3) in the distribution sense. 

Theorem 1.2 (Critical thresholds in 4D case, [10]). Suppose n = A, \ < 0, 
po G D^, and vq G D*"*"-^ with vo{0) = for some positive integer s. Let 
C{r) = VQ{r) + \\\r~'^ pQ{.s).s^ds. The classical solution to (1.4)-(1.7) is 
global if and only if both of the following conditions hold for all R > 0; 

(1) vo{R) ^Oif J^pois)s^ds = 0; 

(2) drC{R) ^ and vo{R) + Rv'q{R) > -y/2RdrC{R); 

where, in the last inequality, we allow the case where the both sides equal 
zero. If Po and vq satisfy the above condition then the corresponding solution 
of (1.4) -(1.7) satisfies 

p€C^{[0,oo),D')nC°°{iO,oo),D'), 

V G C\[0, oo),D'+^) n C"^((0, oo),D'+^), 
$ G C2([0, oo),D'+^) n C~((0, oo),D'+^). 
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The solution is unique in C'^{[0,oo), D^) x C^([0, oo), x C^{[0,oo), D^) 
and also solves (1.1)-(1.3) in the distribution sense. 

Remark 1.3. The above two theorems are true also in the case s = 0. How- 
ever, in that case, p is not differentiable. Then, introducing a new un- 
known m{t,r) := p(t, s)s"^~^ds and replacing (1.4) with the equation 
dtm + vdrtn = 0, we say the solution (p, v, is "classical" in the sense that 
(m,v,$) solves this equation and (1.5)-(1.7) in the classical sense. This is 
also true for Theorems 1.6, 1.7, and 1.12, below. 

Remark 1.4. The assumption fo(0) = is natural because when we try to 
reconstruct the solution (r, v, P) of (1.1)-(1.3), the velocity v should be 

v(t,rc) = {x/\x\)v{t, \x\). 

Remark 1.5. If (p, is a solution to (1.4)-(1.7), then, for any constant 
c, v,^* + c) also solves (1.4)-(1.7). Therefore, in above theorems, the 
solution is unique under certain condition on #, such as <I>(0, 0) = 0. This is 
valid for all results below (Theorems 1.6, 1.7, 1.12 and Corollaries 1.17 and 



One-dimensional case is so simple that everything is made explicit. On 
the other hand, what is special in four-dimensional case is that we can 
write down some integral quantity explicitly. Essentially, their method gives 
necessary and sufficient condition for all other dimensions. However, when 
we try to state them in terms of the slope of the initial velocity, complex 
descriptions are inevitable. This is why they give only sufficient conditions 
for global existence and finite-time breakdown. They mentioned in [10] 
that some further tedious calculations may enable us to obtain a complex 
criterion. 

1.3. Main results. The purpose of this paper is to perform the "further 
tedious calculation" and determine the necessary and sufficient condition for 
global existence/finite-time breakdown in the case other than n = 1,4. We 

introduce a new quantity with which wc state the necessary and sufficient 
condition for global existence (we have already used in Theorem 1.2). For 
n ^ 3, we define 



This quantity represents the balance between the initial velocity and the 
strength of the force governed by the Poisson equation. This quantity clar- 
ifies the conditions for higher dimensions. Note that 



which contains the information about v'q. When we restrict ourselves to the 
case vq > 0, then the use of drC does not change the representation of the 
conditions so much. However, in this paper, we modify the method in [10] 
and allow the case vq ^ 0. The point is that some condition is given in 
terms of C and independent of the sign of vq- For example, if n ^ 3 then 
one sufficient condition for finite-time breakdown is that there exists R> 



1.20). 
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such that drC{R) < (Theorem 1.7, below). Thus, the use of C makes the 
statement slightly clearer. 

Before treating the repulsive case A < 0, wc first illustrate the result in 
the attractive case A > 0. Then, it turns out that the quantity C(r) is very 
useful for the description of the necessary and sufficient condition for the 
global existence. 

Theorem 1.6 (Critical thresholds for attractive case). Suppose A > 0, 

n ^ 1, po E , and vq G D^^^ with Vo{0) = for an integer s ^ 1. 

(1) // n = 1 or 2 then the solution to (1.4) -(1.7) is global if and only 
if po{r) = 0, VQ{r) ^ 0, and drVo{r) ^ holds for all r ^ 0. In 
particular, if pQ ^ then the solution breaks down in finite time. 

(2) // n ^ 3 then the solution is global if and only if 

Voir) ^ 0, C{r) ^ 0, and drC{r) ^ 
hold for all r ^ 0. 

If Pq and vq satisfy the condition for global existence, then the corresponding 
solution of (1.4) -(1.7) satisfies 

p £ C^i[0, oo), D') n C7°^((0, oo), D'), 

V G C\[0, oo), D'+^) n C°°((0, oo), D'+^), 

$ G C\[0, co),D'+^) n C°°((0, oo), D^+2). 

The solution is unique in C'^{[0,oo),D^) x C^{[0,oo),D'^) x C^{[0,oo), D^) 
and also solves (1.1)-(1.3) in the distribution sense. 

In this case, the global existence of the solution is completely character- 
ized as the non-negativity and non-decreasing property of the quantity C. 
Roughly speaking, if C{R) is negative then the attractive force is so strong 
that the characteristic curve starting at r = i2 reaches r = in finite time. 
If drC{R) is negative then it implies that the shock is formed because outer 
wave propagates slower than inner wave does. The proof appears somewhere. 

We now turn to the repulsive case A < 0. In spite the fact that C becomes 
always nonnegative, the situation becomes more complicated; negative vq is 
allowed, that is, the property vq < does not necessarily lead to finite- 
time breakdown. Moreover, positive drC does not necessarily gives the 
global existence, either. We introduce the notion of pointwise condition for 
finite-time breakdown {PCFB, for short) which is a necessary and sufficient 
condition for finite-time blowup given only with the information of initial 
data at r = i?. Rigorous definition is given in Definition 2.8. We also 
introduce the quantity A: 

We now state our main result. 

Theorem 1.7. Suppose X < 0, n ^ 3, pQ e D"^ , and vq G with 
Vo{0) = for an integer s ^ 1. Then, the classical solution of (1.4) -(1.7) 
breaks down in finite time if and only if there exists R such that one of the 
following PCFB (given in Propositions 1.8, 1.9, and 1.10, below) is met. On 
the other hand, the classical solution is global if and only if, for all r > 0, 
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the PCFB does not hold. Moreover, if the condition for global existence is 
satisfied, then the corresponding solution satisfies 

p G C2([0,oo),L>") nC°^((0,oo),L''), 

V e cn[o,oo),L>^+i)nC°°((o,oo),D^+i), 

$ e C2([0, oo),D'+^) n C~((0, oo),D'+^). 

FuHhermore, it is unique in C^{[0, oo), Z>°) x C^([0, oo),D^)x C^{[0, oo),D^) 
and also solves (1.1) -(1.3) in the distribution sense. 

As stated above, this theorem holds also for s = 0, see Remark 1.3. 

Proposition 1.8 (PCFB for vq > 0). Suppose X < 0, n ^ 3, andvo{R) > 0. 
Then, the PCFB is that either one of following three conditions holds: 

(1) drC{R) < 0; 

(2) drC{R,) = and 

drA{R) p 

dy <0; 



vo{R) 2 Jr (C(i?)-^(ii)y-(-2))3/2 
(3) < drC{R) < drA{R)R-(''-^') and 



voiR) 2Jr (C(i?)-A(i?)y-(«-2))3/2 

Proposition 1.9 (PCFB for vq = 0). Suppose A < 0, n ^ 3, arac/ t'o(-^^) = 0. 
Then, the PCFB is that either one of following three conditions holds: 

(1) drC{R) < 0; 

(2) drC{R) = and 

(a) n = 3; 

(b) n = 4 and v'q{R)R < 0; 

(c) n ^ 5 an(f i;(,(i?)i? < - ^''~^^/^ (l - /„), 
where In is a constant given by 



In ■= ((1 - y"')~^ - l) dy 



<i; 



(3) drC{R) > and 
(a) n = 3 anc/ 



VC f RdrC\^ ( VC + VC + RdrC' 

+ ^ 1 T^TT- log 



2c 7 ° \ /sa;^ 



(b) n = 4 anti w^i? ^ -V2RdrC; 
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(c) n ^ 5 and 



{n-2)\{RdrC)l f_, , {n-2)C \^ 



(n - 2)C2 _ RdrC 



2 \ 2C 

Here, we omit R variable in C , drC , and v'^, for simplicity. 

Proposition 1.10 (PCFB for vq < 0). Suppose A < 0, n > 3, and 
vo{R) < 0. Then, the PCFB is that A{R) = or either one of following five 
conditions holds: 

(1) drC{R) < 0; 

(2) drC{R) = a,nd 

' ^ V ;y dy<2drt^{R)- 



\vo{R)\ 2 J J, (C(i?)-^(i?)y-(-2))3/2 
(3) < drC{R) ^ drA{R)R-^''-^^ and 



\vom 2Jr {C{R)-A{R)y-i^-^)fn 
(4) < drC{R) < drA{R){R-^''-^) +VQ{Rf/A{R)) and 

1 1 f^a^cm^ drC{R.)-OrA{R)y-^''--^ ^ , , 



\v^{R)\ 2} J, {C{R)-A{R)y-(--^)f/^ 

(5) +^;o(i?)Vyl(i?)) ^ drC{R), 

where 

'^[cm) dz 



These conditions are very complex but explicit. Once we know the initial 
density po and the initial velocity vq, wc can calculate the condition. Of 
course, in the four-dimensional case, the condition obtained by Propositions 
1.8, 1.9, and 1.10 becomes the same one as in Theorem 1.2. 

Corollary 1.11. If n = A, the PCFB given in Propositions 1.8, 1.9, and 
1.10 is reduced to the following condition: 

(1) A{R) = andvQ{R) < 0. 

(2) drC{R) < 0; 

(3) drC{R) = and vo{R) + v'q{R)R < 0; 

(4) drCln,) > and voiR) + v'q{R)R ^ -^j2RdrC{R). 

In the two-dimensional case, the quantities A and C have different defi- 
nitions. We introduce 

A{r) := 2|A|mo(r) 



and 

C(r) := vo{r)'^ — A{r) log r. 
With these quantities, we obtain a similar theorem. 

Theorem 1.12. Suppose X < 0, n = 2, and po G D^, and vq G D*"*"^ with 
vo{0) = for an integer s ^ 1. Then, the classical solution of (1.4) -(1.7) 
breaks down in finite time if and only if there exists R such that one of 
the following PCFB (given in Propositions 1.13, 1.14, o,nd 1.15, below) is 
met. On the other hand, the classical solution is global if and only if, for 
all r > 0, the PCFB does not hold. If the condition for global existence is 
satisfied, then the corresponding solution satisfies 

p € C2([0, oo),D-') n C7°^((0, oo),D-'), 

V £ C\[0, oo), D'+^) n C°°((0, oo), D'+^), 

$ G C\[0,oo),D'+^)nC°°{{0,oo),D'+^). 

Furthermore, it is unique in C^([0, oo), D°) xC^([0, oo), D^) xC^([0, oo), D^) 
and also solves (1.1) -(1.3) in the distribution sense. 

Proposition 1.13 (PCFB for vq > 0). Suppose X<0, n = 2, and vq{R) > 
0. Then, the PCFB is that 

A{R) 



v'oiR) < 



2Rvo{R) 



eyip{—drC{R)/drA{R)) > R) and that either one of following conditions 
holds: 

(1) po{R) = (drA{R) = Q); 

(2) drA{R) > and 



( arC(R)\ 

1 r'^Pl-OTJ) j drC{R)+drA{R) logy 

MR)^^Jr {C{R)+A{R)logyr/''' 



Proposition 1.14 (PCFB for vq = 0). Suppose X < 0, n = 2, and vq{R) = 
0. Then, the PCFB is A{R) > and that either one of following conditions 
holds: 

(1) po{R) = (drA{R) = 0); 

(2) drA{R) > and 

„ JA(R)RdrA(R) MR) 

RVoiR) ^ - ^ 2 ^ RSrAiR) 

MR) 

2A{R) - RdrA{R) r^^^^^ dz 

4 7i VWz' 

Proposition 1.15 (PCFB for vq < 0). Suppose X < 0, n = 2, and vo{R) < 
0. Then, the PCFB is A{R) = or either one of following conditions holds 
(we omit all R variables, for simplicity): 

(1) PO = (drA = 0); 

(2) drA > and 

(a) dr{v^) ^ A/R + {drA)log{Re^/^); 



(b) A/R ^ dr{vl) < A/R+{drA)\og{Re^'^) and 



1 , 1 r'^K-m drC + drA\ogy rao^,^ 

ay ^ max(U, zarfc*); 



\vo\ 2Jr (C + Alogy)3/2 

(c) dr{v^) < A/R and 



1 /•e^p(-i^) drC + drAlogy 
R {C + Alogyf/^ 



where 



i: 



1.4. Some applications. 

Example 1.16. In the following cases, (1.4)-(1.7) has a unique global solu- 
tion, and the solution solves (1.1)-(1.3) in the distribution sense. 
(1) n = 1, A < 0, and 



Po{r) = e ^ Voir) = \/ , sinr. 



(2) n = 2, A < 0, and 



/'o(^) = T~> — 2' ^o('') = V"-^- 
1 + 

Corollary 1.17. LetX>0orn^3. Suppose po £ D° n ((0, oo), r"~Mr) 
is noi identically zero and ;,'o G satisfies wo(0) = and ^ as r ^ cx). 
T/ien, i/ie solution of (1.4) -(1.7) is global if and only if X > and ra ^ 3, 
and the initial data is of particular form 



/ 2A f 
^°("^ = V (n-2)r-^ i ^o(^)^"-^d.. 

If X > 0, n ^ 3, Po & D'^ n L^((0, oo), r"~^dr) /or s ^ 1, and vo is as above, 
then Vo G D^'^^ and the corresponding solution satisfies 

peC\[0,oo),D')nC'^{{0,oo),D'), 

V e C\[0, oo),D'+'^) n C°°((0, oo),D'+'^), 

# G C2([0, oo), n ^00(^0^ oo), L'^+2). 

Furthermore, it is unique in C^([0, oo), D°) xC-'^([0, oo),D^)x C^([0, oo), D^) 
anti also solves (1.1) -(1.3) in the distribution sense. 

Proof. In the case where n = 1, 2 and A > 0, we deduce from Theorem 
1.6 that the solution breaks down in finite time because po is nontrivial. 
Let n ^ 3, then the assumptions po G L^{{0, oo),r^^^dr) and ^ as 
r — > oo imply C{R} as R ^ oo. Since C(0) = by assumption, we 
see from Theorems 1.6 and 1.7 that the solution is global only if C = 0. In 
the case A < 0, C = implies po = 0, which is excluded by assumption. In 
the case A > 0, the solution is global if we take the positive root vo{R) = 



(n-2)R' 
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Remark 1.18. In Corollary 1.17, the global solution has an explicit represen- 
tation. This is because we can solve (2.4), below, explicitly by separation of 
variables. In this corollary, the case A < 0, n = 1 and the case A < 0, n = 2 
are excluded. If A < and n = 2 then it is not clear whether or not the 
assumption of Corollary 1.17 leads to nonexistence of global solution, but 
following another non-cxistcncc result holds. On the other hand, the case 
where A < and n = 1 must be excluded since the first example in Example 
1.16 is a counter example. This example also suggests that the following 
different version also fails if n = 1. 

Corollary 1.19. Let A < and n ^ 2. Suppose po G is not identically 
zero and vq G satisfies vq{Q) = 0. Suppose, in addition, that there exists 

a sequence {rj}j^i with rj oo as j ^ oo such that VQ{rj) = for all 
j ^ 1, limsnpj^^ rjVQ{rj) < oo, and r'jPo{rj) as j —>■ oo. Then, the 
solution of (1.4) -(1.7) breaks down infinite time. 

Proof. In the n ^ 3 case, VQ{rj) = leads to 

2|A|(rXr,) - (n - 2) po{s)s-^ds) 



drC{rj) = 



(n - 2)r"- 



Since po is nontrivial, J^^ po{s)s'^~^ds > for large j. Moreover, r"po('"j) 
as j — *■ oo by assumption. Hence, we conclude that drC{rj) < for large 
j, which is a sufficient condition for finite-time breakdown. 

Let us proceed to the two dimensional case. We now show that, if j is 
sufficiently large, then the PCFB for R = rj (given in Proposition 1.14) is 
satisfied and so the solution breaks down in finite time. Since po is nontrivial, 
we can suppose A{rj) = 2|A| J^^ po{s)sds > 0. The case po{rj) = is trivial 
and so we now suppose drA{rj) > 0. It suffices to prove that the inequality 



(1.13) rAin) i - ^iSSMSI,^ 



A(rj) 

^2A{rj)-rjdrA{rj) f^^^^^ dz 



\/logi 



is true for some j. Since the left hand side is upper bounded for large j, 
by assumption, it suffices to show that the right hand side is arbitrarily 
large for large j. Notice that the right hand side of (1.13) can be written as 

{Airj)/2)fiA{rj)/rjdrAo{rj)), where 



;(.)= 1 ^+r_^_^r 

■'^ ' A 2xA 



dz 



\/log^ 



Since /(1/2) = -V2iand/' (.t) = {2x^)-^ jf {log z)-^/^dz ^ (2a;5/2)-i(e^- 
1), we see that f{x) ^ oo as x ^ oo. By assumption, A{rj)/rjdrA{rj) = 
Jo^ Po{s)sds/rjpo{rj) ^ oo as j — oo. Thus, the right hand side of (1.13) 
goes to infinity as j ^ oo. □ 

Corollary 1.20. Suppose n \, pq = pc > ^ is a constant, and vq = 0. 
Then, the solution of (1.4)-(1.7) is global if and only if X < 0. //A < then 
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the corresponding solution satisfies 



p G c2([o,oo),L'°°) nc°°((o,oo),i:'°^), 

V G C\[0, oo), L>°°) n C°°((0, oo), D°^), 
$ € C2([0,oo),L'°°)nC°^((0,oo),L'~). 

FuHhermore, it is unique in C'^{[0, oo), L>0) x Ci([0, oo), D^)xC^{[0,oo), D'^) 
and also solves (1.1) -(1.3) in the distribution sense. 

Proof. We first consider positive A case. Since po is not zero, solution breaks 
down if n = 1, 2. In the case n ^ 3, we have C{R) < for all R> 0, which 
immediately leads to finite time breakdown. 

Let us show that the solution is global if A < 0. The one-dimensional 
case is obvious from Theorem 1.1. In the two-dimensional case, we apply 
Theorem 1.12. The PCFB is given by Proposition 1.14 for all -R > because 
vo = 0. Notice that A{R) = \\\pcr'^ > and drA{R) = 2\X\pcr > for all 
R > 0. Therefore, in the end, wc sec that the solution breaks down if and 
only if there exists i?o > such that 

MRq) 

^ 2A{Ro)-RodrA{Ro) /-e WP^ST 

However, the left hand side is zero, and the second term of the right hand 
side is also zero by the relation 2A{R) — RdrA{R) = 0. Since the first term 
in the right side is negative, such Rq does not exist and so the solution to 
(1.4)-(1.7) is global. 

We proceed to the case n ^ 3. The proof is the same as in two-dimensional 
case. Notice that drC{R) = 4\X\pcR/n{n — 2) > and so that Proposition 
1.9 gives the PCFB. In the case n = 4, it is obvious that there does not exist 
i?o such that Vq{Ro)Ro ^ — \/ 2RQdrC {Rq) . In the cases n = 3 and n ^ 5, 
by using the fact that C{R) = 2\X\pcR^/n{n - 2) and so RdrC /2C = 1, we 
verify nonexistence of Rq for which the PCFB holds. □ 

The rest of paper is organized as follows: We first collect some preliminary 
results and illustrate the strategy for proof in Section 2. The main issues 
there are a reduction of the Euler-Poisson and an introduction of the notion 
of pointwise condition for finite-time breakdown. Then, we prove our main 
theorems in section 3. 



2. Preliminaries 

2.1. Reduction of the Euler-Poisson equations to an ODE for char- 
acteristic curve. We reduce the above system (1.4)-(1.7) by employing the 
characteristic curve X defined by an ODE 

^X{t,R) =v{t,X{t,R)), X{0,R)=R 
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and introducing the "mass" 

pr 

m{t,r) := / p{t, s)s"-~^ds. 
Jo 

Then, an integration of (1.4) yields 

(1.4') dtm + vdrm = 0, 

which is written as 

(2.1) jmit,Xit,R))=0. 
Integrating (1.6) and combining with (1.5), we also have 

(2.2) ^X{t,R) - -v{t,X{t,R)) - - • 

Note that (2.1) implies that the mass is conserved along the characteristic 
curve. Thus, we get an ODE for X: 

(2.3) X"(t,i?) = -^^^^^, riO,R)=voiR), XiO,R) = R, 

where mo is the "initial mass" mo(-R) = pois)s"'^^ds. This reduction is 
the key for our analysis. Multiply both sides by X' to obtain 

(2.4) {X'it, R)f = voiRf - , L + '^"""^^^ 



(n-2)i?"-2 {n-2)X(t,RY-'^ 
if n ^ 3 and 

(2.5) {X'it, R)f = voiRf + 2Xmo{R) log 

if n = 2. 

We now state the result about existence of X. We regard X{t,R) as a 
function x — ^ M. For a nonnegative integer s, we define 



C([0,oo)) ifs = 0, 

C([0,oo)) nC^((0,cx))) ifs>0. 



For nonnegative integers si, S2 and intervals Ii, I2, wc define 

C'^'^^ih X h) = {fit, x):hxl2-> md^d'j e C{h X h), 
VaG [0,si],V6€ [0,S2]}. 

Proposition 2.1 (Existence and regularity of solution of ODE (2.3)). Sup- 
pose n ^ 1 and A G M. Let s be a nonnegative integer and assume po G 
and vq G with fo(0) = 0. Then, mo G holds, and for any R > 

there exists t{R) > such that X{t,R) is defined from ODE (2.3) in an 
interval [0,t{R)). Moreover, if there exists T > such that X{t,R) > 
holds for all {t, R) G [0, T) x (0, 00), then we have 

X G c2'^+^([o,r) X (o,oo))nc~'^+i((o,r) x (0,00)). 
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2.2. Local existence of the solution of (1.4)— (1.7). We introduce the 
indicator function 



The interpretation of r(t, R) wih be clear from the following lemma. 

Lemma 2.2 (Lemma 5.1 in [10]). Consider the Euler-Poisson equations 
(1.4)-(1.7). Let X he characteristic curve, then 



Using the above representation of the solution, we deduce the following 
Proposition, which plays a crucial role in our proof. 

Proposition 2.3 (Corollary 5.2 in [10]). The smooth solution to the radial 

Euler-Poisson equations (1.4)-(1.7) is global if and only ifT{t,R) is positive 
for all t ^ and R ^ 0. On the other hand, the smooth solution to the 
Euler-Poisson equations breaks down at t = tc if and only if the following 
equivalent conditions are met for some R = Re'- 



Even if it is possible to determine a function X which solves the ODE 

(2.3) for large time, we can define the solution to the Euler-Poisson equations 

(1.4) -(1.7) by Lemma 2.2 as long as X and V = OrX are positive. 

Proposition 2.4 (Local existence of the solution of (1.4)-(1.7)). Suppose 
n ^ 1 and A G M. Let s be a nonnegative integer and assume po € 
and vq € D^^^ with vq{0) = 0. Let X be the solution of (2.3) given by 
Proposition 2.1. Define T by (2.6). If X{t,R) > and T{t,R) > hold for 
all R> and t G [0,r) and if lim mi R^o^t, R) > for t e [0,T), then 
X{t,0) = for t G [0,r) and (1.4)-(1.7) has a unique solution 



Remark 2.5. In above proposition, if s = then p is not spatially differ- 
entiable. In that case, we use the mass m instead of p and consider the 
modified equations (1.4') and (1.5)-(1.7). 



(2.6) 





(1) /o*=9,t;(r,X(r,i?e))dT 

(2) r{tc, Re) = 0; 

(3) dRX{te,Rc) = Q. 



= — oo; 



pG c2([o,r),L>^)nc°°((o,r),L»*), 

v G c^([o,r),L>^+i)nc°°((o,r),L>^+i), 
$ G c2([o,r),L>"+2)nc°°((o,r),L>^+2)_ 
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Remark 2.6. If (p,v,^) is a solution to (1.4)-(1.7), then, for any constant 
c, (p, $ + c) also solves (1.4)-(1.7). Therefore, in above theorems, the 
solution is unique under certain condition on such as $(0,0) = (See, 
also Remark 1.5). 

Proof. We first show X{t, 0) = for t G [0, T). Since liminfR^o r(t, R) > 0, 
we have R ^ CX(t, R) for small R. Then, the fact that mo{R) = 0{R'^) as 
R^ gives 

I w/.. _ |A|mo(fi) moCR) _ r,(m 

as i? ^ 0. Taking the limit ^ 0, we obtain X"{t,0) = for t G [0,r). 
By X'(0,0) = t;o(0) = and X{0,0) = 0, we have X'{t,0) = X{t,0) = 
for t G [0, r). It gives the continuities of X, X' , and X" (and higher time 
derivatives) around R = 0: 

X G C2'°([0, T) X [0, oo)) n C'^''+\[0, T) X (0, oo)) 

nc°°'°((o,T) X [o,oo))nc°°'^+^((o,r) x (o,oo)). 

Then, the existence part is an immediate consequence of Lemma 2.2. 

Wc prove the uniqueness. It suffices to show in the case s = 0. Let 
{pi,Vi,^i) {i = 1,2) be two solutions to (1.4'), (1.5)-(1.7) which satisfy 

Pi€C\[0,T),D^), 

Vi€C\[0,T),D'), 

^ieC\[0,T),D^). 

Without loss of generality, we can suppose that ^i{0, 0) = since, otherwise, 
we only have to replace ^i{t, x) by ^i{t, x) = ^i{t, x) — $j(0, 0). Now, solving 
■^Xi(t,R) = Vi{t, X(t, R)), we can define the characteristic curves Xi and 
X2, and the indicator functions Fi and Then, we have 

Xi G C\[0,T),D'), Ti G C2([0,r),C((0,oo))). 

Since two solutions exist until t < T, for all i? > and 6 > there exist 
positive constants ci = ci {R, S) and C2 = C2 {R, S) such that 

Xi{t, R)^ci>0 and Tj (i, R) ^ C2 > 0, Vi e [0,T - S]. 

Recall that both Xi and X2 solve 

X''(t,i?) = -A ^|^°^^j_, , X'{0,R) = vo{R), X{0,R)=R. 

Wc fix i? > and 5 > 0. If n = 1 then wc immediately obtain Xi{t, R) = 
X2{t, R) for i G [0,T — 6]. Let us proceed to the case n ^ 2. Using the fact 
that 

77 — 1 

^^^\Xiit,R)-X2it,R)\ 



Xi{t,R)^-^ X2{t,R)^- 
for all t G [0,T — S], and applying Gronwall's lemma to 

Xit, R) = R + X\T)dT, X\t, R) = voiR) - j^'^-i dr, 
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we deduce that X[(t,R) = X2{t,R) and Xi{t,R) = X2{t,R) hold for t G 
[0,r — S]. Since i? > is arbitrary, we also have Xi{t,0) = X2{t,0) for all 
t € [0,T — 5] by continuity. Thus, we see that Xi{t,R) = X2{t,R) for all 
R^ and t G [0, T) since > is also arbitrary. Applying Lemma 2.2, we 
conclude that pi = P2, vi = V2, and so $i = <1>2. □ 

2.3. Pointwise condition for finite-time breakdown. Let us proceed to 

the discussion on global existence. By means of Lemma 2.2 and Proposition 

2.3, it is clear that the existence of tc > such that T{tc, R) = implies the 
finite-time breakdown of the solution. The next elementary lemma suggests 
that the existence of tc > such that X(tc,Rc) =0 with some Rc > also 
leads to the same situation. 

Lemma 2.7. Let X be a characteristic curve. If X{tQ, Ri) = X{tQ,R2) for 
some to > and ^ i?i < R2, then there exist t G [0,to] '^''T'd R G [i?i,i?2] 
such that T{t,R) = 0. In particular, if X(to, Rq) = for some to > and 
Rq > 0, then there exist t ^to and R ^ Rq such that T(t, R) = 0. 

By Proposition 2.3, to ensure the existence of the global regular solution, 
it suffices to start with the initial data for which 

X{t, R)>0, Vi? > and r{t, R) > 0, Vi? ^ 

hold for all t > 0. Now, we introduce the notion of pointwise condition 
for finite-time breakdown, which is "almost" the necessary and sufficient 
condition for the existence of tc G (0, 00) such that r{tc, R) = 0: 

Definition 2.8. For fixed R > 0, we call a necessary and sufficient condition 
for the existence of tc G (0, 00) such that X{tc,R) = or T(tc,R) = hold 
as a pointwise condition for finite-time breakdown. In the case of R = 0, we 
regard a necessary and sufficient condition for the existence of tc G (0, 00) 
such that r{tc, 0) = as a pointwise condition for finite-time breakdown. 
We denote PCFB, for short. 

With this notion. Propositions 2.3 and 2.4 are reduced as follows: 

Proposition 2.9. The local solution to the radial Euler-Poisson equations 
(1.4) -(1.7) given in Proposition 2.4 breaks down in finite time if and only if 
there exist some R^ such that the PCFB is met. 

The meaning of "pointwise" will be clear with this proposition. Therefore, 
specification of the PCFBs is the key for the proof of Theorems 1.7 and 1.12. 

2.4. Construction of the solution of (1.1)— (1.3). At the end of this sec- 
tion, we confirm that the solution of (1.4)-(1.7) solves the original equation. 

Proposition 2.10. Suppose n ^ 1 and A G M. Assume po £ '^f^d 
vq G with vo{0) = 0. Let {p,v,^) be a solution to (1.4)-(1.7) given 
in Proposition 2.4. Then, r{t,x) := p{t,\x\), v(t,x) = (x/|x|)f (t, |x|), and 
P{t,x) := $(t, |x|) solve the Euler-Poisson equations 

(1.1') rt + div(rv) = 0, 

(1.2') vt + V • Vv = -AVP, 

(1.3') AP = r 

in the distribution sense. 
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Proof. Suppose that the solution of (1.4)-(1.7) exists for t < T. We take a 
test function (p{t,x) G C^([0,r) x M"). Since (m,u,$) solves (1.4'), (1.5)- 
(1.7) in the classical sense, the triplet (r,v,P) solves the (l.l')-(1.3') in 
the distribution sense in {\x\ 7^ 0}. Indeed, for all R > e > and fixed 

< t < T, we have 

/ {vt + div(rv))(pdx = / dw dr{mt + rdrm){r)(p{r,Lo)dr 

Je<\x\<R JS'^-'^ Je 

[{mt + rdrm){R)(p{R, to) — {mt + rdrm){e)(p{£, co)] dco 

Sn-l 

— dco {mt + rdr'm){r)drip{r,u!)dr = 0. 

Hence, we only consider the case where ([0, T) x {0}) H supp^^ 7^ 0. Put a 
positive small number £ > and set Q{e) := {x G M"] supj^^j^„ |a;i| < e) 
and Xe{x) := 1q(£)(x). Then, denoting // fgdxdt by {f,g), we have 

(2.10) {rt + div(rv), V?) = {xsrt, ^) + (x£div(rv), 9?) . 

We show that the left hand side is equal to zero. The first term of the right 
hand side is bounded by 

C I r{t, x)dx ]dt + C / r(0, x)dx. 

J \JQis) J JQie) 

We write the inverse map of i? 1-^ X{t,R) by X R{t,X). This is well- 
defined as a map from R+ to itself and the two limits i? — ^ and X — 
are equivalent since dRX{t, R) > for all -R > and X{t, 0) = as long as 
solution exists. By the formula of p given in Lemma 2.2, it holds that 

/ r{t,x)dx ^ / r{t,x)dx 

JQie) J\x\4:V^e 

^ Rit,xr-^pom,x)) ^n-i, 

Jo x^-^r{t,R{t,x)) 

= / poir^-^dr^O 
Jo 

as e — > since po{0) < 00 by continuity. The second term in the right hand 
side of (2.10) is written as J2i=i {Xedi{r\),(p). We estimate only the case 

1 = 1 and n = 2: 

{Xediirv),(p) = J {{rv(p){t,£,X2) - irv(p){t,-£,X2))dx2dt 

/ {rvdi(p){t,x)dxdt. 

JQis) 

Since v{t,0) = ■^X{t,0) = and v is continuous, we see that v{t,x) is 
bounded in [0,T] x Q(e), and so is v. Thus, both terms of the right hand 
side of (2.10) tend to zero because Jq^^^ r(i, x)dx as £ 0, as checked 

above. Therefore, the left hand side of (2.10) is zero since £ is arbitrary. 
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'dx 



Similarly, we can verify that (v^ + v • Vv + AVP, ip) = and (AP — r, (^) = 
0. We only note that r""^^^^ = Jq ps"'~^ds by (1.6), and so that 



/ \VP\dx [ ps^-^ds 

J\x\<e Jo 



as £ ^ 0. □ 

3. Proof of Theorems 

3.1. Proof of Theorems 1.1 and 1.2. Let us first introduce the proof of 
critical thresholds in n = 1, 4 by Engelberg, Liu, and Tadmor in [10]. What 
is special in these cases is that the equation (2.3) can be solved explicitly. 
We begin with the one-dimcnsional case. 



Proof of Theorem 1.1. Integrating (2.3) twice, we immediately obtain 
and so 



V{t,R) = l + v',{R)t+\-^^^t\ 

The solution is global if and only if these two values stay positive for all 
positive time. X{t,R) > holds for all t > if and only if vq{R) ^ or 
vo{Rf - \X\Rmo{R)/2 < 0, and X{t, R) > holds for alH > if and only 
if v'o{R) > or (v(,(i?))2 - \X\po{R)/2 < 0. Therefore, the solution is global 
if and only if 

vo{R) > -^y2\X\Rmo{R), and v'o{R) > -^2\X\po{R) 

for all R> 0. Moreover, it is easy to check that the case vo{R) = mo{R) = 
and the case v'o{R) = p{R) = is also admissible. □ 

Let us proceed to the four-dimensional case. 

Proof of Theorem 1.2. Plugging (2.3) to (2.4), we see that 

(X'it, R)f = voiRf - + = C{R) - Xit, R)X"(t, R), 

which implies {X{t,RY)" = 2C{R). Then, integrating twice gives 
X{t, R) = + 2vo{R)Rt + C{R)t^ 

and so 

^ 2R + 2{vo{R) + v'o{R)R)t + drC{R)t^ 
^' ' 2X{t,R) 

Since ivo{R)Rf - C{R)R^ = -\X\mo{R), X{t,R) > holds for ah i > if 
mo{R) > or if mo{R) = and vo{R) ^ 0. On the other hand, r{t, R) stays 
positive for all positive time if and only if one of the following conditions 
holds 

(1) drCiR) ^ and vo{R) + Uo(i?)i? ^ 0; 

(2) drCiR) > 0, vo{R)+v'q{R)R < 0, and {vo{R)+v'o{R)R)'^-2RdrC{R) < 
0. 

Therefore, we obtain the stated criterion. □ 
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3.2. PCFBs for repulsive n ^ 3 case. We prove Theorem 1.7. From 
Propositions 2.9 and 2.10, all our task is to determine the PCFB, that is, to 
show Propositions 1.8, 1.9, and 1.10. 

Proof of Proposition 1.8. We first note that, by (2.3) and the assumption 

A < 0, X"{t, R)>0 holds as long as X{t, R) > 0. Since X'{0, R) = vo{R) > 
0, we have X'{t,R) > 0, at least for small time t € [0,To]. Note that 
X'{t,R) > for t G [0,To] implies that, for t £ [0,To], X{t,R) ^ X{0,R) = 
R > and so X"(t, R) > 0. Then, it means that X' is also increasing for 
t G [0, To]. Thus, we can choose Tq arbitrarily large, that is, X'{t,R) > 
for all t ^ 0. Then, for all t ^ 0, it follows from (2.4) that 



R \/C{R) - A{R)y-i^-'^) ~ 



This identity tells us that X{t, R) ^ oo as t ^ oo (This also follows from 
the fact that X'{t,R) ^ X'{0,R) = voiR) > 0). For simplicity, we omit the 
R variable in the following. Differentiate with respect to R to obtain 



r{t) 1 1 /•^(*) drC - 



vo 2 7^ 



dy = 0. 



VC - AX(t)-("-2) vo 2Jr [C - Ay-i^-"^))^'^ 
We put 

B{t) := ^ ' = = — + - — "-^ —dy. 

Assume drC{R) < 0. Then, since X(t) ^ oo as t — oo, 

d , drC -drAXit)-^''-^^ ,,,, , drC 
dt 2(C-AX(t)-("-2))'/2 2C3/2 

holds for sufficiently large t. Hence, we have B{t) — oo as t ^ oo, and 
so there always exists a time tg ^ such that i?(to) ^ 0. We see that 
drC{R) < is a sufficient condition for finite-time breakdown. 

Next we assume drC{R) = 0. Then, B{t) is monotone decreasing because 

—B(t) = — :n^X'(t) ^ 0. 

dt ^' 2(C-^X(t)-("-2))3/2 

Therefore, there exists a time io ^ such that B{to) ^ if and only if 
limt_,.oo B{t) < (including the case limj^oo B{t) = — oo). This condition is 
equivalent to 

1 1 /■°° 9,.^!/- ("-2) 



o / o/n-dy<0. 

VQ (C- Ay-("-2))2/2 

Wc finally assume drC{R) > 0. We first consider the case (^|^)"^ > R. 
Then, B{t) takes it minimum at a time t = ti ^ such that 
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because ^B{t) is as above and ti is the time such that ^B{ti) = 0. There- 
fore, there exists a time to such that Bito) ^ if and only if 



1 

1 . 

monotone increasing. Therefore, B ^ -B(O) = I/uq > for all t ^ 0. □ 



We finally consider the case (l^)''-^ ^ R- However, in this case, B is 



Remark 3.1. The argument for above proof is essentially the same as in [10]. 
However, this argument is not directly applicable to the case vq = 0. This 

is because the differentiation of {C — Ay~^^^'^^)^^/'^dy produces the term 
I/vq. Therefore, more delicate analysis is required if ?;o = 0. 

Proof of Proposition 1.9. First note that we have, at least in a small time 
interval, X{t,R) > because X(0,i?) = > 0. Since X"{t,R) > holds 
as long as X{t,R) > by (2.3), we can find a time to > such that 
X'{to,R) > X'{0,R) = vo{R) = 0. Note that to can be chosen arbitrarily 
small. Then, repeating the argument as in the proof of Proposition 1.8, we 
see that, X'{t,R) ^ X'{to,R) > for all t ^ to, which shows X'{t,R) > 
for alH > and X{t,R) ^ oo as t ^ oo. Moreover, X{t,R) ~ C{R)^fH 
for sufficiently large t since X'{t,R) C{R)^^'^ as t — >■ oo. It reveals that 
if drC{R) < then the characteristic curves must cross and so the solution 
breaks down in finite time by Lemma 2.7. 

We now suppose drC{R) ^ 0. We omit R variable in the following. Since 
X'{t) ^ for all t ^ 0, an integration of (2.4) gives 



By a change of variable z = y/R, the left hand side is equal to 
"^W/^ Rdz 



i: 



Vc-^ii:-("-2)z-("-2)' 



We temporally assume that Vq > and take the limit fo i later. This 
computation is justified, for example, by replacing vq by X'{eR,R) > 
with small e > and taking the limit £ J, 0. Differentiation with respect R 
yields 

Q_ RdRiX{t)/R) ^ /-^W/^ dz 



i: 



VC- AX(t)-("-2) Jl Ai?-(«-2)z-("-2) 



^ x{t}/R _ i^drAR-^''-'^) - (n - 2)^i?-("-i))2-('*-2)^^ 



2(C-^i?-("-2)z-("-2))^/^ 



For simplicity, we omit t variable in X and BrX for a while because the 
following computations do not include any differentiation. An elementary 
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calculation shows 
(3.1) 

Q_ QrX _ X ^ /-^/^ dz 



i: 



2C ii (C- Ai?-('»-2)z-{n-2))3/2 



2C 



1 (C- Ai?-("-2)z-{n-2))3/2 



dz 



2(C-Ai?-("-2)z-("-2))3/2 



- AX-("-2) i?VC - AX-("-2) il VC- Ai?-("-2)z-(»-2) 

2C ii ^^_^^_(„_2)^-{n-2))V2 

(-a^CAi? + C5^Ai? - (n - 2)AC) f^/^ ^-("-2) 

—I— / H 

2Ci?-2 A (^_^^_(„_2)^_(„_2))3/2 - 



It also holds that 

-drCAR + CdrAR - (n - 2) AC _ / v'^A drAR - (n - 2) A 

2Ci?"-2 ~ \CR^ ^ 2Ci?"-2 ' 

Now, let us show that 

rX/R ^ 
(^•2) ^0 ((7_^i?-(n-2)^-(n-2))3/2^^ = Ai?-("-2) (n - 2) ' 

Fix a small e > 0. Then, we have 

lim Wo / rrrciz = 0, 

v^iQ ^^_^^_(„_2)^-(„-2))3/2 

since the integral is uniformly bounded with respect to vq. Moreover, 

/•!+£ ^-(n-2) 

■"0 / 77T<^2 

A (C- Ai2-("-2)^-{n-2))3/2 

^ Ai?-("-2)(n-2) ii 2(C-Aii-(«-2)z-("-2))3/2 
< 2^;o(l + e) 



Ai?-("-2)(n-2) 
2(1 + £) 



(C - ~^ -(C- (1 + £)-("-2)) ■ 



1 - 

" 2 



Ai?-("-2)(n-2) 

21 



as vq — > 0. Similarly. 

"l+e _2-(«-2) 

(C-Ai?-("-2)z-("-2))^/^ 



^0 J ; — ; — r-^"^ 



^ 2^0 /■^+^ Ai?-("-2)(n-2)z-("-i) 

^ Ai?-("-2)(n-2) A 2(C-^i?-("-2)z-(-2))='/2 ^ 



^ii;-("-2)(n-2) 

as tiQ — 0. Since £ > is arbitrary, we obtain (3.2). 
Taking the limit i in (3.1), 



CV2yi - {RjXY--^ RC^^^l - (i?/X)"-2 A CV2yi-^-(n-2) 



2C3/2 A ^1 _ ^-(n-2) (n - 2)C ■ 

Thus, we have 

dRX{t) X{t) /-^W/-^ dz 



VI - (i?/X(i))"-2 i?^! - (^A(*))"-' Vl - 2-(^-2) 



2C Ji Vi _ z-{n-2) (n-2)CV2' 

We denote this by B{t). 

Case 1. We first assume that drC{R) = 0. Wc put 

s dz 



G{s) := , " - 

/-. _ („_2) 



An elementary calculation shows, for s > 1, 

(n _ 2'ls~("~2) 

«'W = - 2(l-.U)3/^ <°' 

and so G is monotone decreasing. Moreover, considering the inverse map of 
z {1 — 2;~("""2))-i/2^ .^^g J^avg 

r dz S-1 /-^ A , 

Jl y/l _ 2;-(n-2) _ s-(n-2) 7(i_^_(„-2))-^ V / 

Therefore, 

1 /""^ / 1 \ 

G{s) = - ({1- y~')-^^ - l) dy. 

One verifies that if n = 3 then limg^oo G{s) = — oo. We now put, for n ^ 4, 

For any m > i ^ 4 and y € (1, oo), it holds that 

(l_y-2)-7^ < (l_y-2)-7^. 
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This gives Im < Ii ^oic m > I A. If n = 4 then 

/■TV 

h= lim / 

N^oo Ji \ / 

= hm ((A^^ - 1)5 - (AT- 1)) = 1. 

Thus, we obtain 

— oo if n = 3, 

hm G{s) = {0 if n = 4, 

1 - /„ > if n ^ 5. 



Since 



B{t) = G{X{t)/R)- '^'"'"^ 



{n - 2)CV2 ' 

we conclude that there exists to £ [0, oo) such that T(to) ^ if and only if 

(1) n = 3; 

(2) n = 4 and v'qR < 0; 

(3) n ^ 5 and v'^R < - /„). 

Case 2. We assume that drC{R) > 0. We write B{t) = H{X{t)/R). 
Then, it holds that 



ds 2(1 - S-(»-2))3/2 2C(l-S-("-2))V2- 

Therefore, the minimum of H, hence of S, is 



RdrC 

The solution breaks down in finite time if and only if this value is less than 
or equal to zero. This leads us to the condition 



(n - 2)C^ 2("-2) 
^ RdrC 



f-, , (n-2)C\ 1 

n-2i RdrC \ rV+^id^) dr. 

-C2 —— 1 ' 



2 \ 2C J Ji ^1 _ ^-(n-2) 

Using the identity 



y/l _ ^-(n-2) U i?a^C J I \ {n-2)C 
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we obtain the equivalent condition 



1 3 . . 

{n - 2)2 {RdrC)2 ( (n-2)C7\2(^ 



AC V RdrC 



(n - 2)C5 RdrC 



2 V 2C 

1 



(n-2)C J 

In particular, if n = 3 or 4, then the above integral is computable, and we 
have more explicit condition 



v',R ^ -l./cTm:c + ^ (1 - ^) log 



'Vc_jWc^±^^^ 



if n = 3 and 

v'oR ^ -^/2RdrC 

if n = 4. □ 

Proof of Proposition 1.10. We first note that if A{R) = 0, then X'{t,R) = 
vo{R) < 0. Therefore, the solution breaks down no latter than t = R/\vo{R)\ 
by Lemma 2.7. Hence, we assume A(R) > 0. Then, since X'(0,R) = 
voiR) < 0, X'{t, R) = -VC- AX(t)-("-2) as long as X'{t, R) ^ 0. Take 



r 



dy 



(A)^ ^C - Ay-("-2) 



1 



c ) dz 



We see that, for aU t G [Q,U), X{t,R) > X{U,R) = (A(i?)/C(i?))^/("-2) > 
and X'{t,R) < X'{U,R) = 0. Since X"{t^,R) > by (2.3), using the 
same argument as in the proof of Proposition 1.9, we have X'{t, i?) ^ for 
all t^U and so 

x'(t R) - l- VciR)-MR)x{t,R)-^'^- ^), fovt^U, 

^' ' \ \/C{R) - A{R)X{t, i?)-("-2) , ioit^U. 

We also obtain X(t, R) ^ 00 ast ^ 00. In the following, we omit R variable. 
For sufficient large t, X{t) ~ C^/^t holds since X'{t) C^/^ as t ^ cx). It 
implies that if drC{R) < then the characteristic curves must cross and 
so the solution breaks down in finite time by Lemma 2.7. Differentiation of 
X{U,R) = (^/C)V("-2) with respect to R gives 



drUX'{U,R) + dRX{U,R) = dr (^^^ 
Using the fact that X'{t*) = 0, we obtain 



A\ "-2 



dRX{U,R) = dr(-^ 



1 

A \ "-2 
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Hence, if dr^A/C)^/^'^ ^-^ ^ then the solution breaks down no latter than 
• 

Thus, we assume drC{R) ^ and {dr{A/C){R)Y/^'^-'^^ > in the fol- 
lowing. Notice that the latter condition is equivalent to the following two 
conditions: 



Step 1. We determine the condition that solution can be extended to 

h< 

^ dy 



time t = t^. For t ^t^, we have 



/ 



t. 



lx(t) yJC - Ay-("-2) 
Differentiation with respect to R yields 



For ^ i < i*, it holds that 

< sjc - AX{t)-^^-'^) ^ \/C - Ai?-("-2) = |t;o|. 
Therefore, r(t) has the same sign as 



dy = Q. 



B (f) = r(t) ^ J_ _ 1 a,c-a,^y-("-2) 

VC-AX-(n-2) \vo\ 2y^(,) (C-Ay-(-2))3/2 



Taking time derivative, one verifies that Bi takes it minimum at t = ti G 
[0,i*) such that 

drA\ "-2 \ 



X{ti,R) = min i?. 



Note that (^/C)^/^'*-^) < by assumption, and that {drA/drCy/^''-^^ < 

R is equivalent to drC > drAR''^'^~^\ Since we have already known that 
r(0) = 1 > 0, the solution can be extended to the time t = U UNLESS 
drC > drAR-^-^-^) and 

\drC ) 

is satisfied. Notice that this condition is a sufficient condition for finite-time 

breakdown. 

Step 2. We consider the condition that the solution can be extended 
from the time f = to t = oo. For simplicity, we suppose that solutions are 

extended to time t = (we keep assuming ^ drC < drA{R~^^~'^^ +Vq/A) 
holds). Recall that, for t ^ U, X'{t) = - AX(t)-("-2) 5> 0. As in the 
case vq = 0, this inequality with X"(t) > gives X(t) ~ C^^'^t ^ oo as 
t ^ oo. 

We define t** as the time that t** > U and = R. Then, we have 

t - l"" 
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Therefore, i** = 2t^ and 



x{t) 



dy 



for all t^tif. As in the prevfous step, we set 



B2{t) :-- 



1 /■^(*) drC - drAy-^''-'^'^ 



dy—2drt^. 



- AX{t)-i^-^) \vo\ 2Jii (C - Ay-in-2))3/2 

B2{t) and T{t) has the same sign for t ^ tt,. We also note that B2{t) oo 
ast iU because T{Q > and ^/C^^AX(f)-^^^ ^ as t [ U. It holds 
that 

d drC-0,.AX{ir("-^ 
dt^'^^^ - 2(C-AX(t)-("-2))3/2^ 

(1) If drC{R) = then is monotone decreasing because ^B2{t) ^ 0. 
Therefore, solution can be extended to t = 00 if and only if 



dy — 2drt^ ^ 0. 



(2) If drC{R) > then B2 takes it minimum at t = ^2 such that X{t2) = 
(9r^/5rC)V("-2). Therefore, solution can be extended to i = 00 if 
and only if 



1 I r\ 

\Vo\ ^ JR 



( drA \ \ 
\drC J 



drC — drAy' 



-(n-2) 



(C-Ay-("-2))3/2 



dy — 2drt* > 0. 



□ 



Before proceeding to the two-dimensional case, let us sec that Theorem 
1.7 gives the same criterion as in Theorem 1.2 if n = 4. Namely, we prove 
Corollary 1.11. 

Proof of Corollary 1.11. Before the proof, we prepare some elementary com- 
putations. We note that 



I Br A 



drC dr-C — drAy 



R {C-Ay-'^)^ 

I V 8rC y 

JR 



dy 



drC 

~c' 

drC 



- Ay/ 



AdrC - CdrA 

Tj2'^y + - 



c 



L 



/ drA 
8rC 



R {Cy^-Afl'^ 



dy 



C^-^-A 



{CR^-A)2 



+ 



AdrC - CdrA 



C2 



{CR^-A)-2 - C 



drC 



A 



'^'^^^-'iCdrA-A8rC)h-\y^drC- ^^^A-AB^C 



C2 



C2 



C'^\vq\R 



26 



and that 



1 _ \vo\R _ CdrA - AdrC 

1 , vlR^ + A^ ^ drA \ \vo\R. 



\vo\^ 2C^\vo\r'^'^ 2C\vo\r) ' 
2CR + R^drC - drA _ \vo\R 

. f V0 + Rv'q vqR \ 



(signi;o)5, 



VqR 



C 

where we have used VqI^ + A = CR? and 
2CR + R^drC - drA 

= {2vlR + ^) + {2v,v',R- + drA - ^) - drA 

= 2vqR{vq + Rv'q). 
It also holds that 



\fAC^ 



A C 

Prom Propositions 1.8, 1.9, and 1.10, we see that drC < is the sufficient 
condition for blow-up. Moreover, the PCFB in the case drC = is 

<i^(c*A - Aa^c)i + * (^) = 5i±j?i < 

if vq > 0, 

Rv'o < 

if vq = 0, and 

^^^^^\cdrA - AdrC)l -dr('^\+ 2dr ( = "° < 



(j2 y r r , \ C J \ C J C 

if vq < 0. Hence, the PCFB is summarized as vq + Rv'q < 0. 

Let us proceed to the case drC > 0. If > then Proposition 1.8 implies 
that the PCFB is drC < drAR'^ ^vq + Rv'q < C/vq and 

(3.3) ^^^{CdrA - AdrC)-^ + dr ^ 0. 

We put a = vq + Rv'q, (3 = voRdrC/C > 0, and 7 = drC{CdrA- AdrC)/C^. 
Note that, by assumption, we have < A/C < R^ < drA/drC, which 
implies 7 > 0. Then, (3.3) can be written as a ^ /? — ^7- We make 
this condition clearer. An elementary computation shows that 5 := 7 + 
2aP -p^ = 2RdrC > 0, and that S - 2aP = R9rCi-drC+drAR-^} ^ ^ 
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latter one means < 7. Thus, the inequahty a ^ P — ^77 < is reduced 
to a ^ -V^7 + 2a/3-/32 = -V^, that is, t;o + ^ -^/2RdrC. This 
condition is stronger than drC < drAR~^ <^ vq + Rv'q < C/vq. 

If drC > and vq = 0, then it immediately follows from Proposition 1.9 
that a ^ -VS is the PCFB. 

We next consider the case drC > and vq < 0. Proposition 1.10 gives 
the PCFB. If drC ^ drAR'^, then the condition is 

^^'^^\cdrA - AdrC)h -dr('^\+ 2dr f 0. 



^2 .^^r^^ ^^^r^j \ C J ' ' \ C ^ 

We keep the above notations a, /3, 7, and 5. Then, this is written as a ^ 
j3 — ^/7. Note that the right hand side is negative. By the same argument as 
above, it is also written as a ^ -V^. IfdrAR-'^ < drC ^ drA{R-'^ + V^/A), 
then the condition is 



^I^{CdrA-AdrC)-2 ^ 



dr 



vqR 
C 



which is written as ^ ^ |a — Note that drC < drA{R~''^ + Vq/A) = 
CdrA/A is equivalent to 7 > 0. By assumption, we also have P < and 
7 — /?2 = (5 — 2aP < 0. We now show that a ^ (3 leads to the contradiction. 
In this case, .y/7 ^ |q: — /?| is equivalent to a ^ (3 + y^. However, this is also 
written as 

0<V7^\a-P\=a-(3 ^ 7 + 2a(3 - = 6 > 

a ^ Vd or a ^ —VS. 

The last inequalities cannot be iquivalent to a ^ /3 + ■^/7 since > and 
P + y/^ < 0. This is the contradiction. Hence, P ^ a. Then, ^/7 ^ |q! — /9| = 
P — a corresponds to a ^ —y/6. 

We finally treat the case drC ^ drA{R'''^+VQ/A). We prove this condition 
is stronger than a ^ —VS. An elementary computation show that drC ^ 
drA{R-^ + v^/A) implies 

C voRdrA 

a ^ \ tTa — < 0- 

VQ 2A 

Moreover, introducing the function P{t) = drCt^ + 2at + 2R, we see that 

^ = ^ K-^) (-^) ^ (-^) 



1 

C2 



/ , drA 2A\ 2„2 

- 2{v^ + Rv'q)voR (^vi + A j + 2i? (^vl + 
2vov'oA + J 



1 



v^drA 
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□ 



3.3. PCFBs for repulsive 2D case. We finally prove the two-dimensional 
case. Though we can calculate the characteristic curve in an implicit way 
([10]), we use the argument similar to the previous n ^ 3 case. 

Proof of Proposition 1.13. We first note that X'{t,R) > vo{R) > 0, Vi ^ 
follows from the same argument as in the proof of Proposition 1.8. Then, 
X{t, i?) — oo as t — oo, and, by (2.5), 

rX{t,R) 



L 



R ^vo{RY + A{R)\og{y/R) 



t. 



for all t^O. For simplicity, we omit the R variable in the following. Differ- 
entiate this with respect to R to get 

m _ 1 _ 1 f""^'^ 2vov'o-A/R + drA\og{y/R) 
X'{t) vo 2 J J, {vl + A\og{y/R)f/^ ^ ' 

We put 

Wt) •= = 1 + i r^*^ 2vov',-A/R + drA\og{y/R) 

X'{t) vo^2Jn {v', + A\og{y/R)r/^ 

Since X'{t) > for all t ^ 0, B{t) and T(t) has the same sign. Since dj-A ^ 
by definition, the right hand side is positive for all time if 2voVq — A/R ^ 0. 
Now, we suppose 2vqVq — A/R < 0. Recall that X{t) — > oo as i — > oo and 
that A and vq are independent of time. If drA = then one sees that there 
exist to > such that 



L 



X{to) 



dy 



R (i;2 + ^log(y/i?))3/2 vo\2vov'^-A/R\ 

since J^^^\vq + Alog{y/R))~^/'^dy — oo as t — oo. This implies T{to) = 
B{tQ) = 0, which lead to finite-time breakdown. Let us proceed to the case 
drA > 0. An elementary computation shows that the minimum of B is 

B I e J , Therefore, under the assumption vqVq — A/R < and drA > 0, 



there exists a time to such that r(to) ^ if and only if 

1,1 r^P(-i^) drC + drAlogy 

vo 2Jr (C + ^logy)3/2^^^"- 



□ 



Proof of Proposition I.I4. Let us begin with pointing out that the exactly 
same argument as in the proof of Proposition 1.9 shows X'{t,R) > for all 
t>0 and X{t, R) ^ 00 as t ^ 00. We omit R variable in the following. As 
in the proof of Proposition 1.9, we temporarily suppose that vq > and let 
■uo — later. Integration of (2.5) gives 



t. 



f 

Jr ^yv^ + Alog{y/R) 
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By a change of variable z = y/R^ the left hand side is equal to 

^^(*)/« Rdz 



i: 



^Jvl + A\og■, 
Hence, differentiation with respect R yields 







RdR{X{t)/R) 



^vl + A\og{X{t)/R) 

-R 



+ 



X{t)/R 



dz 



^/^^'+A]ogz 
-dz. 



2{v^ + Alogzf^^ 
For a while, we omit also t variable. An elementary calculation shows 
(3.4) 







X 



^v'^+A\og{X/R) R^v^ + Alog{X/R) 
RdrA f^l^ dz RvldrA r^/^ 

Ji M+Aiogz ^ y 1 



+ 



dz 



y/vf+A\ogz 
dz 



2A ./i ^v^ + Alogz 
dz 



+ Alogz) 



3/2 



RvqV\ 



r 



+ Alog z) 



3/2 ■ 



We now show that 



(3.5) 



lim vq 



X/R 



dz 



{vl + A\ogz) 



Fix a small e > 0. Then, we have 

cX/R 

lim Vq 



dz 



3/2 A' 



0, 



,3/2 



voiO Ji+e {yl + A\ogzY 
since the integral is uniformly bounded with respect to vq. Moreover, 

"1+^ dz 



Vq 



1 {vl + A log 

2vQ{l + e) 
A 

2^o(l + £) 
A 



3/2 



A 



2z{vl + ^logz)3/2 



dz 



Vq 



- {vl + A\og{l + e))- 



2(1 + £) 
A 



as vq — 0. Similarly, 

"1+^ dz 



Vq 



i: 



{v^ + Alogz) 



3/2 



l+e 



A 



^VQ 

A J I 2z{v^ + Alogzf/^""^ ' A 



dz 



as Vq — 0. It proves (3.5) since £ > is arbitrary. 
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Taking the limit i in (3.4), 

Q_ dnX X/R ^ 1 [^/^ dz 



Thus, we have 

dRX{t) X{t)/R /-^W/-^ dz 



B{t) :- 



2A Ji AV2- 
Case 1. We first assume that drA = 0. We put 

s dz 



y^logiX{t)/R) ^logiX{t)/R) Ji 



G{s) := 



-i: 



An elementary calculation shows G'{s) = — (l/2)(log s)"^/^ < for s > 1, 
and so G is monotone decreasing. We also see that G' is not integrable, and 
so that lims_>oo G{s) = —00. Since 

B{t) = G{X{t)/R) + 

we conclude that there always exists to £ (0, 00) such that T{to) = 0. 

Case 2. We next assume that drA > 0. We write B{t) =: H{X{t)/R). 
Then, it holds that 

d , 1 RdrA 
-tH{s) = -— + 



ds ' ' 2(logs)3/2 2A(logs)V2- 

A 

Therefore, the mininiTim oi H , hence of B, is H{e^<>rAy xhe solution breaks 
down in finite time if and only if this value is less than or equal to zero. This 
leads to the condition 

VnR ^ - J eR9rA + — I dx. 

2 V 2VAJJ0 

□ 

Proof of Proposition 1.15. If A{R) = 0, then X'{t,R) = vo{R) < for all 
t ^ 0. Therefore, we deduce from Lemma 2.7 that the solution breaks down 
no latter than t = R/\vo{R)\. Hence, we assume A{R) > 0. Then, since 
X'{0,R) = vo{R) < 0, X\t) = -^yvoiR)^ + A{R)log{X{t,R)/R) as long 
as X'{t,R) ^ 0. Put 

^ dy R r""'" dz 



> re 



Re-o/^ y/v^ + A\og{y/R) ylV2e^oM J I Vlop' 

Then, one sees that, for t G [0,U), X{t,R) > X{U,R) = Re'^'o/^ > and 
X'{t, R) < X'{t*,R) = 0. Since X"{U,R) > 0, the same argument as in the 
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proof of Proposition 1.9 shows that X'{t, R) ^ for all t ^ t^.. and so that 



^ [ - VvoiR)^ + A{R) log(X(t, R)/R ), for t ^ 
^' ^ [ V^o(i?)^ + A{R) log(X(t, fort^f,. 

X{t, i?) ^ oo as i — > oo is also deduced. We omit R variable in the following. 
Differentiation of the identity X{t^,R) = Re~'"o/^ with respect to R gives 

a«X(t.) = e-^oM(l-i?a,(|)). 

Hence, if Rdriv^/A) ^ 1 then the solution breaks down no latter than 

t = t*. Thus, we assume Rdr{vQ/A) < 1 in the following. This is equivalent 

to drv^ < A/R + {v^/A)drA and to -C/A < -drC/drA. 

Step 1. We first consider the condition that solution can be extended to 
time t = t*. For t ^t^, we have 

dy 



I 

Jx 



= t. 



Ix{t) \/C + A\ogy 
Differentiation with respect to R yields 

1 V{t) 1 drC + drAlogy 



--/ 



dy = 0. 



VC + A\ogR ^C + A\ogX{t) 2 Jx^t) {C + A log 2/)3/2 
For ^ i < U, 

< ^yC + A\ogX{t) ^ ^yC + A\ogR = \vo\ 
holds. Therefore, r{t) has the same sign as 

B (f) = r(t) ^ ^ _ 1 drC + drAlogy 

y/C + A\ogX{t) \vo\ 2 7^(,)(C + ^logy)3/2 

Taking time derivative, one verifies that Bi takes it minimum at t = ti G 
[0,U) such that 

/ ^ drC^ 

X{ti) = min I R, exp 



drA^ 

Here, note that = exp{—C/A) < X{ti) by assumption. Also note 

that ^ 

exp ( — TT— 7 ) ~ ^exp 



f drC\ „ / drV^-A/R\ 



Since we have already known that r(0) = 1 > 0, the solution can be extended 
to the time t = U UNLESS drV^ > A/R and 

„ , 1 1 drC + drA]ogy 

is satisfied. Notice that this condition is a sufficient condition for finite-time 

breakdown. 

Step 2. We next consider the condition that the solution can be extended 
from the time t = to t = oo. For simplicity, we suppose that solutions 
are extended to time t = U (we keep assuming drV^ < A/R+ {vQ/A)drA 
holds). Recall that, for t>U, X'{t) = y/C + A\ogX{t) > 0. 
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We define t^^ as a time such that i** > and X{T^^) = R. Then, we 
have 



and so = 2t^. Thus, 



Ir ^C + Alogy 
for all f^t*. As in the previous step, we set 



t 



r(t) 1 , 1 /-^w 5,c + 5,^iogy , 



r(t) 1 1 r 

JC + A\ogX(t) \vo\ 2 7^ 



^/C + A\ogX{t) \vo\ ^Jr (C + ^logy)3/2 

_ 1 1 /-^W drv',-{A/R)+drA\og{y/R) 
~\vo\^2Jn (C + Alogy)3/2 

i?2(i) and r(t) has the same sign for t>U. Wc also note that i?2 ^ oo as 
i i because r(i*) > and y/C + A log X ^ as t J, i*. It holds that 

d drvl-{A/R)+drA\og{X{t)/R) , 
Jt^'^^' = {C + AlogX{t))V^ ^ 

If dj-A{R) = then B2 is monotone decreasing by assumption drV^ — A/R < 
0. Moreover, ^B2{t) is uniformly bounded by (9r^o ~ (^/-^))/l^o| < from 
above, and so there exists time t2 such that .62 (t2) = 0. Therefore, now we 
suppose drA{R) > 0. 

B2 takes it minimum at i = i2 such that X{t2) = exp{—drC/drA). There- 
fore, the solution can be extended to t = 00 if and only if 

„ . 11 r^P(-S^) drC + drAlogy , ^ . 

1^0 1 2Jr {C + Alogyy/^ 

□ 
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